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Abstract 

The peculiarities of electric current in semiconductors with nonuniform distri- 
bution of charge carriers are studied. The semiclassical drift-diffusion equations 
consisting of the continuity equations and the Poisson equation are solved nu- 
merically using Rusanov finite-difference scheme of third order. The different 
types of boundary conditions are numerically investigated. It is shown that the 
stability of the Rusanov scheme for the problem considered is achieved with the 
Neumann type boundary conditions. These conditions correspond to the absence 
of diffusion through semiconductor surface. Special set of parameters is found 
under which a very interesting and unusual transient effect of negative current 
in nonuniform semiconductors appears. Different regimes of negative current are 
considered for realistic semiconductor materials. 
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The study of electric current in semiconductors is important for describing and 
modelling semiconductor devices. The action of external radiation can result in the 
formation of a nonuniform distribution of charge carriers in a semiconductor sample. 
The transport processes in semiconductors with essentially nonuniform distribution of 
carriers can be quite specific. For instance, in a semiconductor sample biased with an 
external constant voltage the electric current may turn against the applied voltage [1,2]. 
Of course, it can happen only as a short time fluctuation after which the current turns 
back becoming positive. So we can talk about this only as about a transient effect. 
This unusual transient effect of negative current is interesting by itself although there 
might be some useful applications suggested. For example, it can be used for measur- 
ing physical characteristics of semiconductor materials or for constructing electronic 
controlling devices [3]. 

Transport properties of semiconductors are usually described by the semiclassical 
drift-diffusion equations [4]. In dimensionless units we have the continuity equation 

where 0<x<l,£>0,i = l,2, and the Poisson equation 

dE 

— = 47r( Pl + p 2 ) . (2) 

Here, p 1 = pi(x,t) and p 2 = p 2 {x, t) are charge densities; E = E(x,t), electric field. 
The functions p 1 , p 2 , E are unknown functions to be found from equations (1) and 
(2). 

The considered semiconductor sample of length L is characterized by the following 
parameters: pi, mobilities; Di, diffusion coefficients; 7*, relaxation widths. Two types, 
(i — 1, 2), of charge carriers are assumed: positive and negative, so that p\ > 0, p± > 0, 
and p 2 < 0, P2 < 0. We consider the semiconductor sample of the length L — 1. 

Initial distribution of charge carriers is modelled by the Gaussian form, which e.g. 
corresponds to the distribution of charges in a semiconductor irradiated by an ion 
beam. For t — 0, 

ft(M) = §J ex p{-^f^} > (° ( 3 ) 

where Qi = Jq pi(x, 0) dx and Zi = Jq exp j — j dx. Note that the initial distri- 
bution for the electric field E(x, 0) is to be found from equations (1) and (2), using (3) 
and given boundary conditions. 
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There are two possibilities of posing boundary conditions for functions Pi(x, t). The 
Dirichlet type, 

p t (0,t) = Pl (l,t) = 0, (^ = 1,2), 

which corresponds to the absence of the drift current (ji = PiPiE = 0) at the semicon- 
ductor surface, and the Neumann type, 



dpi 
dt 



dpi 



dt 

x=0 



0, (^ = 1,2), 



which corresponds to the absence of diffusion through the semiconductor surface. 

We considered both types of boundary conditions for Pi . It turned out that calcula- 
tional procedure is less stable in the case of the Dirichlet boundary conditions because 
of the accumulation of charges at the semiconductor surfaces. As a result, discon- 
tinuities appear in the vicinity of the boundaries x = and x — 1, which leads to 
numerical instability. In practice, it is also more difficult to arrange that the drift cur- 
rent through a semiconductor surface be absent. In what follows, we use the Neumann 
type of boundary conditions for the system (1) and (2). 

The role of boundary conditions for the electric field E(x, t) is played by the voltage 
integral 

/ E(x,t)dx = l. 
Jo 

Our aim is to study the peculiarity in the behaviour of the electric current through 
the semiconductor sample 

J(t) = I jtot(x,t) dx , (4) 
Jo 

where 

(d \ ( d \ 1 dE 

^-D 1 -y 1 +^E-D,-y + --. ( 5) 

We also are interested in the electric currents at the left (x = 0) and right (x — 1) 
surfaces of the semiconductor, J (t) = jt o t(0, t) and J x (t) = j tot (l,t). Let us note that 
J(t) can be presented in the integral form 

J(t) = [ \uipi(x,t) +fx 2 P2(x,t)] E(x,t) dx+ 
Jo 

+D 1 [pi(0, t) - Pl (l, *)] + D 2 {p 2 (0, t) - p 2 (l, *)] , 
which is more convenient for numerical computation. 
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Numerical investigation showed that the known explicit finite difference schemes 
of the second order are strongly unstable for the system (1) and (2). Because of this, 
we resorted to the third-order accuracy finite-difference scheme suggested by Rusanov 
[5] which proved to be very efficient for discontinuous solutions. The calculational 
procedure for the Rusanov scheme is constructed as follows. We consider the vector- 
function vffl = {pi(x m ,t n ), p 2 (x m ,t n ), E m (x m ,t n )}, where x m = mh and t n = nr, 

m = 0, 1, ... , M; n — 0, 1, The grid space step is h — 1/M. The time step r is 

chosen so that the calculational procedure be stable. The parameter a in the Rusanov 
scheme, responsible for the stability of calculations, is a = 0.45 for r/h = 0.125. We 
will not present the computational scheme explicitly since it is quite cumbersome. Let 
us better dwell upon the approximation of boundary conditions for pi(x,t) (i = 1,2) 
and for E(x,t). Numerical modelling showed that the Neumann type of boundary 
conditions for the charge carriers are to be approximated with the third order accuracy 
for the time layers t — t n — nr and auxiliary time layers t = tff — t n + 2r/3 at the 
integer boundary nodes: 

pi{x M ,t) = ^ [4pi(x M -i,t) - Pi(x M -2,t)} , Pi(x ,t) = ^ [4p i (x 1 ,t) - pi(x 2 ,t)} , 

and with the second order accuracy for the half-integer boundary nodes at the auxiliary 
time layer t — t$ = nr + r/3: 

Pi(xM+l/2,t) = Pi(x M -l/2,t) , Pi(X-i/2,t) = Pi(x 1/2 ,t) , 

where x m +i/2 = mh + 1/2, m — 0,1, ... , M, n = 1, 2, . . .. 

To approximate the boundary conditions for the electric field E(x,t), given in the 
integral form, it was necessary to rewrite it in the form convenient for calculations: 



E(x,t) = 1+4:71 



rl rx 

Q(x,t)- Q(x,t) dx , Q( x ,t)= [pi(y,t) + p 2 (y,t)] dy 
Jo J Jo 



Then E(0, t) and E(l, t) are given by the above formulas for x = and x — 1. At each 
time layer t n = nr, we checked whether the equality Jq 1 E(x,t)dx = 1 holds true. 

We have accomplished numerical calculations for different sets of parameters a^, 6j, 
Pi, Di, 7j, and Qi. It turned out that for some values of parameters there exists an 
unusual transient effect when the electric current turns against the applied voltage. 
This effect of negative current occurs on the whole manifold of parameters. Here we 
illustrate it for the case of a unipolar semiconductor, when Q 1 = Q > and Q 2 = 0. 
Other parameters are set as p\ — 1, p 2 — 0, Dj <S 1, 7, 1, and the notation a\ = a, 
b\ = b is used. 
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Fig. 1 shows the total electric current (4) through the semiconductor sample, as a 
function of time, for the injected electric charge Q = 0.5, the initial Gaussian distribu- 
tion width b = 0.1, and for different locations of the peak of this charge distribution, 
a. As is seen, the negative electric current happens at the initial stage of the process 
for a < 0.5, such as a = 0.1 and a = 0.25, but does not exist for a > 0.5. 

Fig. 2 presents the total electric current versus time for Q = 3 and b = 0.5, with 
varying locations a. Again, the negative current appears for a < 0.5 but does not occur 
for a > 0.5. More accurately, the conditions on a and Q, sufficient for the occurrence 
of the negative current, are 

11 ^ 1 

a< 2-^Q> Q> 2^- 

The existence of the transient effect of negative electric current is the most inter- 
esting result of the numerical modelling we have accomplished for semiconductors. 
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Figure captions 



Fig. 1. Total electric current as a function of time for Q = 0.5, b = 0.1 and 
different locations for the peak of the initial charge-carriers distribution: a = 0.1 (1 
curve, solid line); a = 0.25 (2 curve, long-dashed line); a = 0.5 (3 curve, short-dashed 
line); and a = 0.75 (4 curve, dotted line). 

Fig. 2. Electric current through the semiconductor sample versus time for Q = 3, 
b = 0.5, and varying initial charge locations a, which take the same sequence of values 
as in Fig. 1. 
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